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^P \ A novel method, based on superpotentials is proposed for obtaining the quasi-normal 

modes of anti-deSitter black holes. This is inspired by the case of the three dimensional 
BTZ black hole, where the quasi-normal modes can be obtained exactly and are propor- 

qh, tional to the surface gravity. Using this approach, the quasi-normal modes of the five 

dimensional Schwarzschild anti-deSitter black hole are computed numerically. The modes 

W)' again seem to be proportional to the surface gravity for very small and very large black 

holes. They reflect the well-known instability of small black holes in anti-deSitter space. 
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1 Introduction 



Numerical studies of perturbations of black holes have revealed the presence of certain charac- 
teristic modes that govern the time evolution of the initial perturbation. These modes seem to 
depend only on the black hole parameters and not on the nature of the perturbation. This was 
first recognised by Vishveshvara [0 while studying perturbations of Schwarzschild black holes. 
Since then, perturbations of asymptotically fiat black holes have been analysed and these char- 
acteristic modes(called quasi-normal modes) have been found for some black holes (a detailed 
account can be found in 0, 0). 

Recently, there has been an interest in quasi-normal modes of AdS black holes in light of 
the AdS /C FT correspondence @]. There is some evidence to suggest that an off-equilibrium 
configuration in the bulk AdS space is related to an off equilirium state in the boundary theory. 
In [§], black hole formation by collapse of a thin shell in AdS space is investigated. The shell is 
related by AdS /C FT duality to an off-equilibrium state in the boundary theory which evolves 
towards equilibrium when the shell collapses to form a black hole. 

A quasi-normal mode governs the decay in time of a perturbation of the black hole con- 
figuration in the bulk. It should therefore be related by AdS/CFT duality to the return of 
the boundary Yang-Mills theory to thermal equilibrium. The role of ingoing modes in deter- 
mining the thermalization time scale of the boundary Yang-Mills theory was pointed out in 
1^0. A numerical computation of quasi-normal mode frequencies for AdS black holes in various 
dimensions has been done in [^. 

In this paper, we compute the quasi-normal modes for the ArfS'-Schwarzschild black hole in 
five dimensions using a superpotential approach. This is done for the perturbation of a black 
hole by a minimally coupled scalar field. Time-independent normal modes of a scalar field in 
the background of this black hole have already been analysed P, ^. 

We notice, that for the (2 + 1) — (i BTZ black hole, exact quasi- normal mode solutions can 
be obtained from the Klein-Gordon equation as the black hole potential belongs to a class of 
exactly solvable potentials, derived from a superpotential. 

We propose, based on analogy with the results for the (2 + 1) — d BTZ black hole, and 
on numerical evidence in five dimensions, that the black hole potential can be described by a 
potential series derived from a superpotential. The problem is then exactly solvable and the 
form of the quasi-normal mode wave function is thus known. This is used as an ansatz to obtain 
the quasi-normal modes of the AdS^ black hole. However, the numerical values of the mode 
frequencies do not agree with the earlier calculation in |^. We comment on this discrepancy. 

■^The normalisation measure used here is not appropriate, however, and these modes are not quasi-normal 
modes. 



The organisation of this paper is as follows: The next section is a brief review of quasi- 
normal modes and their properties, and a short summary of some numerical approaches to 
computing these modes. 

The actual computation of quasi-normal(QN) modes for the AdSs-Schwarzschild black hole 
is described in Section 3. First, the QN modes for the (2 + 1) — d BTZ black hole are evaluated. 
Motivated by this, a superpotential approach to finding the QN modes for the five- dimensional 
black hole is presented and the mode frequencies obtained. 

In Section 4, we discuss our results. We show that the behaviour of the modes computed 
numerically by us for small black holes is consistent with that expected from the differential 
equation obeyed by the mode solutions. We mention salient points of two earlier papers P and 
on the subject and comment on the discrepancy between the numerical values obtained by 
us and those in 0]. We also state work in progress on QN modes for the RN AdS^ black hole, 
and black holes in other dimensions. 



2 Quasi-normal modes and their properties 



Quasi-normal modes of a black hole are characteristic modes associated with the decay of any 
perturbation outside a black hole. In general, these modes do not form a complete set, in 
the sense that the time decay of a perturbation cannot be described completely in terms of 
them. However, they dominate the decay at certain intermediate or late times. To compute 
these modes, one studies the decay of an initial perturbation by imposing ingoing boundary 
conditions at the horizon and (for asymptotically flat black holes) outgoing boundary conditions 
at infinity. This ensures that no gravitational wave from the horizon or infinity disturbs the 
initial perturbation. Due to these boundary conditions, the quasi-normal modes are complex. 
The mode wave functions are non-normalisable in space but exponentially decay in time. Since 
the exact mode wave functions are not known even for the Schwarzschild black hole, there exist 
only numerical computations of quasi-normal modes. These are also very difficult to do as it 
is not easy to isolate the purely ingoing wave near the horizon that is also outgoing at infinity. 
Such a wave function blows up at these points and could be contaminated by a small outgoing 
part at the horizon or a small ingoing part at infinity where these parts go exponentially to zero. 
A unique boundary condition to isolate quasi-normal modes has been given by Nollert . The 
QN modes are poles of the Green's function corresponding to the wave equation obeyed by the 
perturbing field. Using this, the QN modes for the Schwarzschild and Reissner-Nordstrom black 
holes were evaluated numerically [[l^], |jll[. WKB methods have also been used to calculate 
these modes [0 (see also [0 for more references). 



The method of continued fractions ||T3| gives very accurate results for the modes. Here, after 



removing the singular parts at the horizon and infinity, the QN mode wave function is expanded 



as a series. The coefficients in the series obey a three-term recurrence relation and depend on the 
quasi-normal mode frequency. The frequencies then have to satisfy a continued fraction relation. 
An interesting analytic approach to computing the QN modes is by approximating the black 
hole potential by an exactly solvable potential (the Poschl Teller or Eckart potential) |jl4| , [|15] . 
Ferrari and Mashhoon [Q have mapped the problem of QNMs of the Poschl- Teller potential 
to the bound states of the inverse potential and computed the modes. These approximate the 
fundamental QN modes of the Schwarzschild black hole well. 

For the case of asymptotically anti-deSitter black holes, the black hole potential diverges at 
infinity. Therefore, as pointed out in |]^, the quasi- normal modes are defined to be those that 
are ingoing at the horizon, but fall off to zero at infinity. It can be shown that this also agrees 
with the unique boundary condition of Nollert for obtaining quasi-normal modes. However, 
due to the divergence of the AdS black hole potential at infinity, most of the methods described 
above to compute the QN modes of asymptotically flat black holes cannot be used directly. It 
is seen that a combination of two methods, inspired by the case of the {2 + 1) — d black hole, 
yields good results. 



3 Numerical computation of quasi-normal modes 



We ffist consider the case of the QN modes of the non-rotating (2 + 1) — d BTZ black hole [|16 
The non-rotating BTZ black hole metric for a spacetime with negative cosmological constant 
A = — p- is given by 

ds^ = - {Nf de + {N)'^ dr^ + r^ {d<pf (1) 

with 

N = \k-M + ^) (2) 



The horizon radius r+ = y Ml. We consider a massless scalar field in the black hole back- 
ground. The Klein-Gordon equation for the scalar field is written using an ansatz (independent 
of the angular variable) for the field 

^ = —^ X{.r) exp{iujt) (3) 

and by going to the tortoise coordinate r^^ where dr^ = -^. The Klein- Gordon equation is 

d'x 



drl 



+ V{r)x = uj'x (4) 



where 



Since r = — r+ coth(ar^,) where 
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a is the surface gravity of the black hole and is equal to 2txT where T is the Hawking temper- 
ature. 



This potential can be obtained from a superpotential [jT^ of the form 

W = A coth(Q;r*) + B tanh(ar*) 



by 



Here, from (0), 



V = W^ - {Wy - A^ - B^ - 2AB 



A = -a{l/2±l] 



B 



a 
2 



(9) 

(10) 
fll) 



As discussed in |[T^, the lowest energy state for the potential W'^ — (W)' has energy zero, so 
the lowest energy for the potential V is —{A + 5)2. The wave function corresponding to this 
state is 



X 



exp(-|(iy)) 



(sinh(ar,))-^/(") (cosh(ar,))-^/(") 



(12) 



For any of the two values of A and B in (pT]), the wave function (|l^) is not normalisable. It 
blows up either at the horizon or at r = oo (i.e r* = 0). We are interested in the quasi-normal 
mode solution that blows up at the horizon and falls off to zero at infinity. This corresponds 
to A = —{3/2)a. Also, B = — ^. The lowest energy (formally, since the wave function 
corresponding to this is non-normalisable) is then 



E = -4(a)2 
and the lowest quasi-normal mode cu = yE. Therefore, 



(13) 



u = -2ia (14) 



Thus, for a given horizon radius r+ = y Ml, we can calculate a taking / = 1. Then we can 
read off the value of the lowest quasi-normal mode. 

The case of the quasi-normal modes of the ArfS's-Schwarzschild black hole is not so simple 
as the black hole potential is more complicated than the BTZ case. 

The metric for this black hole is 

ds^ = - {Nf de + (Ar)-2 dr^ + r^ {dn^f (15) 



with 



N = ^1 + (r2)/(/2) _ (ro/r)2 (16) 

where the black hole mass M is 

and A3 is the area of a unit 3-sphere. The Klein-Gordon equation for the scalar field in the 
background of this black hole can be written, as before, as a potential problem. Using the 
ansatz for the field 

-IN 3/2 

$ = ( - y(r) exD(iujt) (It 



( -j x{r) exp{iujt) 



and by changing to the tortoise coordinate given by dr^^ = ■^, the Klein Gordon equation is 



d'^X , T./„^ _ 2 



drl 



+ Vir)x = oj'x (19) 



where 



^'« - (T-i + l^)(i + ^^-(7'^) (^°) 



Here, for simplicity, we have taken / = 1. We have also not considered an angular dependence 
for the field, which is a simple generalisation of the following analysis. 



If the potential could be written as an exactly solvable potential in the r* coordinate, then 
the quasi-normal modes could be obtained easily from the 'lowest' energy, as in the BTZ case. 
As mentioned in the previous section, the fundamental QN modes of the Schwarzschild black 
hole were closely approximated by the QN modes of the Poschl- Teller potential. This is also 
an exactly solvable potential derivable from a superpotential [|1^ . 



We propose that the potential can be written as a series, derived from a superpotential of 
the form 

oo 

W = 5] A„ coth(nar,) (21) 

n=l 

where as in (0), a is the black hole surface gravity. 



dr ^ ^' ^ ' n 



a = (1/2) ^^^(r+) = r^ + 1 (22) 



As before, a = 2'kT, where T is the Hawking temperature. 
The potential derived from this series as in (^) is 



^1 sinh^(nar^) ^ ^^J^^^.^^ tanh(nar^) tanh(mQ;r*) 

It has lowest energy 

oo 

E = -(E^n)' (24) 

71=1 

The form of the potential is such that it goes exponentially to zero in r,, as r -^ r+ and 
blows up as 1/rl as r* — i> (or r -^ oo). This reproduces the behaviour of the black hole 
potential. We propose that this potential series is exactly equal to the black hole potential 
(pO|). Terminating this potential series at finite order would then mean an approximation to 
the black hole potential in the spirit of the Poschl- Teller method for asymptotically fiat black 
holes. The wave function corresponding to the lowest energy state is 

/oo 
(W)) = n(sinh(nar,))-^"/("°) (25) 

n=l 

It is non-normalisable, as in the (2 + 1) — d case. In order to calculate the quasi-normal modes 
for the AfiS'-Schwarzschild black hole, we use (^) as an ansatz in the equation (|19|). However, 
we truncate the product in the ansatz upto some finite order. This is equivalent to truncating 
the potential series. We first write 

X = i' exp(-2Co'r*) (26) 



This isolates the non-normahsable part of the wave function. We can rewrite (|1^) as an equation 

for ip as 

- ^ + 2tuj^ + Vir)^ = (27) 

drf dr^ 

Here, V is the black hole potential (^). We now try to solve for uj by taking the ansatz 



^ 



■0 = 11(1 ~ exp(2na;r*)) "q (2^ 



ra=l 



This ansatz is obtained from (^) which was an ansatz for x- Thus, A^ represents the order 
upto which Vs has been truncated. 



The ansatz for ifj is substituted into (|27|) and the l.h.s of ( ^Tl) is expanded as a series in 
(x — x+) where 

1 

X = - 
r 

x^ = ^ (29) 



From (p7D, each term of the series is to be equated to zero. We first divide ( |27|) by ip. Then, 



we have 



N N 

- 4(^ AraK^f - 4 a (^ mA^Kra) - 

m=l m=l 

N N 



where 

exp(2mar,) 

^m = 777^ r (31) 

1 — exp((2mar*) 



Expanding each of the terms in the l.h.s of (|30D , 

K„r = Y^^rnin) (x - X+)" (32) 

n=0 

Kl = Y.drn{n) (x-x+)" (33) 



n=0 

V = ^ K (x - x^ 

n=0 



(34) 



Substituting these expansions into (|30|), and equating the term of order (x — x^) in the Lh.s to 



zero, we obtain 



/15 , 3< 9r^. 

Al = - ^\ " 7 . ^' (35) 

2(«w - a)ci(l) ^ ^ 

Similarly, equating the term of order {x — x+)^ in the l.h.s to zero, we obtain a recursion 
relation for Aj^. 



^ {Na - tiu)cNiN) {Na - zcj)c^(A^) 

_ {J2m,n=l Y.N1,N2;N1+N2=N A-mAnCm{Nl)Cn{N2)) 

{Na - iuj)cN{N) 



{Na - iuj)cN{N) A{Na - iuj)cN{N) 



(36) 



The An coefficients are functions of uj and of the black hole parameters. The recursion 
relation for the An is complicated, unlike the case of the method of partial fractions [^ applied 
to asymptotically fiat black holes. There, an ansatz was taken for the wave function (removing 
its singular part) and it was expanded as a series whose coefficients obeyed a simple three-term 



2 



recursion relation. From (27), as r^, — > (i.e as r -^ oo), the two kinds of solutions are ip ~ r* 

5 

and ip ^ ri . The first solution is not normalisable. We want the quasi-normal mode solution to 
vanish as r,, ^ 0. So we choose the second solution. But from our ansatz ( pSj) for the solution, 
this implies a relation between the coefficients An- More precisely, 

^ A 
E - = - 5/2 (37) 



Since the An coefficients are functions of u, this relation gives the value of the lowest quasi- 
normal mode as a function of the black hole parameters. This would then be an approximation 
to the actual QN mode of the black hole at order A^. Thus, the asymptotic behaviour of the 
actual solution is explicitly demanded of the ansatz. This is used to obtain the mode frequencies 
from the relation (131 



We have used this method to calculate the QN modes at various orders N. We find that 
as we increase N, there is a convergence in the mode frequency. The mode frequencies are 



given in Table 1. They have been calculated at order N = 18. Beyond this order, as (|37|) is 

9 



Radius r+ 


Re(tu) 


Im(cj) 


1 


0.6948 


1.4648 


2 


1.0713 


1.9817 


5 


2.4462 


4.2642 


10 


4.8249 


8.3279 


50 


24.0159 


41.3183 


100 


48.0251 


82.6165 


150 


72.0358 


123.9190 


500 


240.1150 


413.0500 


750 


360.1720 


619.5740 


1000 


480.2290 


826.0980 









Table 1. 



complicated, it becomes computationally time-consuming to find its roots. However, as there 
is a clear convergence in the mode frequency as one increases the order from A^ = 1 and as 
the percentage difference between the real and imaginary parts of the mode frequency at order 
iV = 18 and A^ = 17 is about 1% (for Reoj, it is 1.26%, and for /ma;, it is 1.09%), these numbers 
should be a good approximation to the QN modes of the black hole (more discussion on results 
is given in the next section). 



4 Discussion of results 



The lowest quasi-normal mode frequencies for different horizon radii r+ are given in Table 1. 
As mentioned in the last section, there is a convergence in the mode frequency as the order A^ 
is increased. The convergence curve for r^ = 10 is given in Fig.l. The rate of convergence does 
not seem to depend on the value of the horizon radius r_|_. 

The real and imaginary parts of the mode frequencies are plotted as a function of r^ for 
large r+ in Fig. 2a and Fig. 2b. 
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Figure 1. Convergence in Im{uj) as order N is increased, for r_|_ = 10 



It is seen that the both the imaginary and real parts of the mode frequency are proportional 
to r+ for large r_|_. The real and imaginary parts of the mode frequency can also be plotted 
as a function of the surface gravity a. It is seen that for both small and large r+, the real 
and imaginary parts of the frequency are approximately proportional to a. Re{uj) ~ 0.24 a 
and Im{ijj) ~ 0.41 a provides a good fit for the data. In the case of the QN modes where the 
modes could be obtained exactly, we saw that they were proportional to the surface gravity. 
Our numerical results seem to suggest that this may be true even for the five dimensional black 
hole - at least for very small and very large black holes. The surface gravity a = (2r+ + — ). 
Therefore, a is large for very small and very large black holes. We have verified numerically that 
for very small black holes, the mode frequencies are very large, and tu ~ — . This behaviour 



of the mode frequency is expected - as can be seen from the differential equation ( pTf ) on 
changing to the inverse radial coordinate x given by (pOf). Then, the differential equation in 
these coordinates is 



-TqX 



+ x^ + x^ 



dx^ 



+ (2x^ -Ar^x^ + 2iujx'^)-r^ 

ax 



,15 3x2 



+ ^)^ 



(38) 
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Figure 2a. Re{u}) as a function ofr_^ for large r_| 



Now, we scale x as x = g x+ (where x+ = — ). Now, (|38|) is 

/ 26 6, 24, 2\^ y 

{-x.q' - g + x,q + q ^ 



{2xW - ^x\(i - Aq^ + 2iijx+q^^ 



dq^ 
di 
dq 

•4 4 4 4^^ 



+ 



(39) 



Then we see that for very large black holes (small x+ approximation), x+ can be scaled away 
from (|39|) near the horizon in this approximation provided oj = — where C is a constant 
independent of a;+. It can also be checked that near the horizon, there are no solutions to 
the scaled equation with C = 0, except the trivial solution. This shows that u is indeed 
proportional to r+ (i.e — ) for large black holes. 



For very small black holes (i.e in the large x+ approximation), again x+ can be scaled away 
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Figure 2b. Im{u) as a function ofr^ for large r_| 



from (|39|) near the horizon provided uj = D x^ where D is a constant independent of a;+. 
Here too, it can be checked that near the horizon, there are no solutions to ( pQ]) in the large 
x+ approximation with D = (except the trivial solution). This implies that for small 
black holes, cu is proportional to — (i.e to x+). This is also just a reflection of the fact that 
a very small ArfS'-Schwarzschild black hole has negative specific heat and actually resembles a 
Schwarzschild black hole. 



We now examine the previous numerical work on QN modes of AdS black holes in four, five 
and seven dimensions. 



The numerical results obtained in 0] for the five dimensional black hole do not agree with 
our values for the QN modes. There, the solution to (pT]) is expanded as a series around the 
horizon in the inverse radial coordinate x = 1/r. Using this, the l.h.s of ( p7|) is expanded 
as a series around the horizon, and a recursion relation is obtained for the coefficients in the 
expansion of the solution. The coefficients a„ are functions of the black hole parameters and 
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Figure 3a. Re(uj) as a function of a 



the mode frequency u. Therefore, 



(40) 



The mode frequencies are then the roots of the equation obtained by setting the series to zero 
at a; = 0. In actual computation, the series is truncated, and the roots obtained. However, the 
mode frequencies seem to go to zero as r_|_ ^ 0, which is not the behaviour expected from the 
differential equation. 
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Figure 3b. Im(uj) as a function of a 

A numerical computation of QN modes is tricky because of the nature of the boundary 
conditions on the mode solution at both the horizon and at infinity. The mode, which is ingoing 
(and not normalisable) at the horizon, could be contaminated by an outgoing component which 
goes to zero there. At the other boundary, the solution is required to go to zero. In a numerical 
truncation, it could also be contaminated as a; — ;► (i.e as r — ;► oo) by the solution that blows up 
at this end. Therefore, in a numerical computation, it is essential to ensure the correct ingoing 
behaviour near the horizon and the correct asymptotic behaviour of the solution that goes to 
zero as r ^ oo. In our method, similar to the continued fraction method for asymptotically 
flat black holes, we have a specific ansatz for the wave function where the behaviour at both 
boundaries is explicitly present in the form of the ansatz. The form of our ansatz ensures that 
there is no contamination from the outgoing part at the horizon. At the other end as a; — * 
(i.e as r -^ oo), we demand that the solution must fall off as x^ at every order. This ensures 
that there is no contamination from the solution that blows up at this end. 

We have presented a new approach to computing the quasi-normal modes of AdS black 
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holes. The novel feature of this method was an ansatz for the QN mode wave function which 
was derived from a superpotential. This was made in analogy with the case of the three 
dimensional BTZ black hole where the modes can be obtained exactly and the wave function is 
derived from a superpotential. The BTZ QN modes were proportional to the surface gravity. 
The modes obtained by us numerically for the five dimensional black hole are also approximately 
proportional to the surface gravity. More importantly, the modes obtained by us numerically are 
proportional to the inverse of the horizon radius for small black holes, reflecting the well-known 
fact that these black holes are unstable. We have also shown, from some scaling properties of 
the differential equation obeyed by the mode solutions, that this is indeed to be expected for 
small black holes. Work is in progress to compute the QN modes for AdS black holes in four 
and seven dimensions, and also the corresponding Reissner-Nordstrom black holes using this 
approach. Some preliminary work on the five dimensional RN AdS black hole seems to suggest 
that the QN mode increases with the charge of the black hole. 



We would like to thank G.T.Horowitz and V.Hubeny for useful comments about the results 
in our paper. 
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